Dynamical decoupling of a singlet-triplet qubit afflicted by a charge fluctuator 
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The efficiency of dynamical decoupling pulse sequences in removing noise due to a nearby charge fluctuator 
is studied for a singlet-triplet spin qubit. We develop a numerical method to solve the dynamical equations for 
all three components of the Bloch vector under a general pulse protocol, where pulses are applied along an 
arbitrary rotation axis. The qubit is shown to undergo both dephasing and dissipative dynamics, pending on 
its working position. Analytical solutions are found for the limits of weakly and strongly coupled fluctuators, 
shedding light on the distinct dynamics in the different parameter regimes. Scaling of the qubit decay time with 
the number of control pulses is found to follow a power law over a wide range of parameters and qubit bias 
points. 

PACS numbers: 03.67.Lx, 03.67.Pp, 73.21.La, 73.23.Hk 



I. INTRODUCTION 

Considerable attention has been given in recent years to the 
proposal to encode the logical qubit states into spin singlet 
(S), and unpolarized triplet (To) states of two electrons lo- 
calized in a double quantum dot (QD). 1 ~- In a Bloch sphere 
representation of the S - T qubit, \S) = ~ l4-t))/V2 
and \T ) = (|t|) + ||t))/ v / 21ie in the north and south poles, 
and rotations around the z and x axes are performed by the ex- 
change interaction, J, and a magnetic field gradient across the 
two dots, Sh, respectively. Concerted efforts have resulted in 
a series of impressive advances in the initialization, readout, 
control, and coupling of these systems. Rapid electrostatic 
control over the interdot bias provides a highly tunable ex- 
change interaction that enables singlet preparation and single- 
shot readout by utilizing Pauli spin blockade, as well as fast 
rotations around the z axisA- Indeed, aside from their robust- 
ness against uniform nuclear fluctuations, the main advantage 
of S — To qubits over single spin qubits is their amenability to 
fast single-qubit operations. 

To complete single-qubit control, nuclear polarization cy- 
cles, in which bias is swept across the S — T + degeneracy 
point, have been used to exchange spin polarization between 
the electrons and the nuclei, 5,6 and were shown to generate 
different Overhauser (hyperfine) fields in the two QDs.— To- 
gether with J, the resulting field gradient has provided an all- 
electric scheme to perform single qubit rotations around an ar- 
bitrary axis. 7 These nuclear pump cycles were later perfected 
by utilizing the hyperfine coupling in a feedback loop, which 
not only generated a stable nuclear field gradient of Sh = 23 
mT, but also produced a narrowed nuclear state distribution, 
leading to a prolonged dephasing time of T 2 * — 94 ns. 8 Other 
methods to generate local magnetic field gradients were also 
demonstrated, including inhomogeneous Zeeman fields gen- 
erated by on-chip micromagnets.- 

Implementing two-qubit gates in a system of two double- 
dots has proven to be a challenging roadblock, and a first ex- 
perimental demonstration of conditional operation on S — Tq 
qubits was reported only last year. 10 Improved gate perfor- 
mance was recently demonstrated in a work that included a 
complete measurement of the system's density matrix using 
state tomography. 11 While both experiments used capacitive 



coupling to generate a CPHASE gate (that can be transformed 
into a CNOT gate with the addition of two Hadamard gates on 
the target qubit), the most notable new feature in the design 
of the latter experiment is the incorporation of nuclear state 
preparation that generated a stabilized field gradient. Shulman 
et al. then used this gradient (Sh = 5 mT, Zeeman energy of 
0.125/ieV) to apply a spin echo (SE) control pulse along the 
x axis, mitigating charge noise and enabling them to work 
near the singlet avoided crossing, where sizable couplings be- 
tween the two double dots can be obtained. At this bias, where 
J ^> Sh, noise due to nuclear fluctuations is suppressed and 
charge noise plays a significant role. The employed SE pulse 
extended the two qubits' coherence time and allowed for their 
entanglement with a Bell state fidelity of 0.72»i! 

Evidently, dynamical decoupling (DD) pulse sequences 
have been applied to S — Tq qubits prior to the above ex- 
periment with remarkable success, starting with a single pulse 
SE, 4 and following by more advanced control sequences, in- 
cluding CaiT-Purcell-Meiboom-Gill (CPMG), Concatenated 
DD (CDD), and Uhrig DD (UDD) schemes^. The co- 
herence time was extended to a record of more than 200^s 
using a 16 pulse CPMG sequence and pulse optimization 
techniques. 14 In these experiments, dephasing times were 
measured at a negative bias, where exchange is small, and 
decoherence is attributed mainly to the varying Overhauser 
fields. The it pulses used in these sequences were thus ap- 
propriately performed (roughly) along the z axis by pulsing 
the bias near the anticrossing where J ^ Sh. The discussion 
above suggests, however, that when the qubit resides closer 
to the anticrossing (e.g., during two-qubit operations), nu- 
clear fluctuations (causing x rotations) will be suppressed and 
charge-noise-induced fluctuations (causing z rotations) will 
become more prominent, thereby requiring tt x pulses to ef- 
fectively extend the qubit coherence. 

The effectiveness of various decoupling schemes at mit- 
igating nuclear-induced dephasing has been extensively ad- 
dressed theoretically for single-spin qubits 1518 and for two- 
spin qubits under SE. 19 In contrast, their effectiveness in han- 
dling charge noise is less clear. In the current work we address 
this question by extending numerical and analytical stochastic 
methods that were developed in the context of superconduct- 
ing qubits. We shall consider a qubit coupled to a single two- 
level-fluctuator (TLF), treating the latter as a classical source 
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of random telegraph noise. While quantum telegraph noise 
was considered before (see, e.g., Ref. |20|). our classical treat- 
ment of the TLF seems reasonable, given the fact that the area 
surrounding the qubit is likely to be depleted of charge traps. 

The qubit Hamiltonian reads H q = B • er, where B = 
i (Sh, 0, J), and a is the vector of Pauli spin matrices for 
the pseudospin states 5 and To. Unless otherwise noted, 
we take 5h = 0.125/ieV, and the qubit working position, 
ip = arctan((5/i/ J), is determined by the interdot bias, e, 
that controls J. Previous studies were mostly focused on 
the </? = point, where pure dephasing is expected, and the 
performance of the pulse sequence depends only on the TLF 
characteristics^^ Bergli and Faoro considered periodic DD 
(PDD) at the ip = ir/2 point, where both dephasing and dis- 
sipative dynamics take placed Here we extend the work re- 
ported in Ref. |24|, by treating an arbitrary working position 
and a general pulse sequence. 

While we present results pertaining for two-spin qubits in 
gate-defined GaAs double dots, our work is relevant for a va- 
riety of systems, including superconducting qubits and other 
QD materials. In particular, 5 — To qubits were recently im- 
plemented in a Si/SiGe double QD, where dephasing time of 
360 ns was measured. 25 We expect that charge noise will play 
a dominant role in Si, where the hyperfine interaction strength 
is three orders of magnitude smaller, due to reduced coupling 
to- and number of nuclear spins, as compared with GaAs. DD 
has been successfully implemented in other systems such as 
electron spins in irradiated malonic acid single crystals. 26 Fi- 
nally, very recently DD was incorporated with two-qubit gates 
in a hybrid system consisting of an electron spin and a nuclear 
spin in a single nitrogen-vacancy center in diamond.— 

The paper is organized as follows. In Sec. II we present our 
model and discuss the qubit-TLF couplings and their depen- 
dence on various system parameters. In Sec. Ill we detail the 
transfer matrix numerical approach to calculate qubit deco- 
herence under a general pulse sequence. Sec. IV provides an- 
alytical results for the PDD and CPMG cases, shedding light 
on the effectiveness of these sequences at different parameter 
regimes. In Sec. V we use our formulation to estimate qubit 
coherence times for various scenarios, and in Sec. VI we pro- 
vide a short summary of our work, and outline future research 
directions. Appendix A includes formulae for the special case 
of pure dephasing, and Appendix B explains how to use the re- 
sults given in Sec. IV to explicitly find the dynamics of a qubit 
with an initial state along the equator of the Bloch sphere. 



H. QUBIT-TLF COUPLINGS 

A TLF residing near a double dot couples differently to 
the two-electron 5 and To spin states due to their different 
charge distributions. 28 This results in fluctuations in the ex- 
change interaction that lead to qubit decoherence and gate 
errors. 29 ' 30 Various sources can contribute to charge noise in 
lateral gated devices, including donor centers near the gate 
electrodes, switching events in the doping layer, and charge 
traps near quantum point contactsi^ 2 - Charge noise measure- 
ments in GaAs QDs revealed a linear temperature dependence 



characteristic of 1// noise, 33 which was shown theoretically 
to emerge from a TLF ensemble with an exponentially broad 
distribution of switching rates. 34-36 These charge fluctuators 
behave classically and are characterized by switching rates, 
7±, and qubit coupling strength, v. In this picture we can 
write the qubit-TLF interaction as 



Hint = v£(t)a z 



(1) 



where £(t) = ±1 is a classical noise representing a random 
telegraph process, switching between ±1 with rates 7±. 

In a previous work we have developed a multipole expan- 
sion technique to calculate the Coulomb couplings between 
the 5 — Tq qubit orbital states and the fluctuator. 30 Assum- 
ing the TLF is a two-site trap, sufficiently remote from the 
double dot, there is no qubit-TLF tunnel coupling, and the in- 
teraction Hamiltonian includes the terms: —vpaf — v 7 a^aj, 
where ctj (crj) is the qubit (TLF) Pauli operator. We iden- 
tify two possible scenarios for qubit-TLF couplings: (i) 7- 
coupled, where the charge fluctuates between two sites in the 
trap, and (ii) /3-coupled, where the charge jumps in and out 
of the trap. /3-coupled TLFs require a nearby charge reservoir 
(such as the 2DEG layer or quantum point contacts), and are 
expected to be less abundant. In the classical limit considered 
here, the qubit-TLF Hamiltonian reads: 



H = H q + Hi nt = B(t) • a, 



where B(t) is 



B(f) 



K<5fc,o,j-2« 7 e(t)), 



7 — coupling 



^(8h,0,J + vp — vp(i(t)), j3 - coupling 



(2) 



(3) 



It is stressed that in the case of classical TLF considered in 
the current work, the two couplings have the same qualitative 
effect on the qubit dynamics, and here we distinct them only 
for the purpose of evaluating their strength. 

The vp and v 1 couplings were calculated in Ref. [3(3 within a 
Hund-Mulliken orbital model using a multipole expansion to 
quadrupole-quadrupole order. It was shown that the leading 
term in the /3-type (7-type) coupling is dipole-charge (dipole- 
dipole), where left and right entries correspond to the qubit 
and TLF, respectively. 37 Analytical approximations of these 
couplings are found by keeping only a subset of the orbital 
two-electron states: 5(0, 2), 5(1, 1), where indicate the 
number of electrons in each dot. Within this simplified model 
we find the leading terms: 



21 / \ C®*q J (^) 

V t> (£) = ^2J 2 ( £ ) + T 2S] 
ca q a t J 2 (e) 
' ~ £3 2J 2 (e)+T c 2 



(4) 



"TO 



[sin 9t cos 4>t — 3 sin 6 cos <p 
x sin 9 sin Ot cos(0 — 4>t) + cos6>cos6>t] , (5) 



where the left (right) superscript denotes contribution from a 
particular multipole moment of the qubit (TLF): charge (1), 
dipole (2), etc. Here, couplings are normalized to QD con- 
finement energy, hu>o, c = (e 2 / nas) /fouo is the Coulomb 
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to confinement energy ratio, a q (at) is the half interdot (TLF 
intersite) separation, R is the qubit-TLF distance, and tilde 
denotes length normalized to the QD Bohr radius as- Fur- 
thermore, J{e) is the bias-dependent exchange and T c is the 
Coulomb-assisted tunnel coupling between the two dots. The 
angular dependence of vp and v~ ( is specified by four angles 
(0, 0, #t,0t), where the first two define the orientation of 
the qubit-TLF axis, and the last two define the TLF intersite 

•30 

axisias 

Fig.Q]shows the dependence of the qubit-TLF couplings on 
various parameters. Here and throughout the paper we have 
diagonalized the full Hund-Mulliken Hamiltonian, which in- 
cludes all two-electron Coulomb matrix elements. Neverthe- 
less, the analytic formulas in Eqs. (@J, and (0 are instructive 
as they qualitatively capture many of the couplings features. 
We model the double dot with a quartic potential with dot con- 
finement hujQ — 3meV (as ~ 20 nm), and a q — 2.8, and use 
the dielectric constant for GaAs, k = 13.1. In addition, un- 
less otherwise noted, we take the TLF center radius D t = 5 
nm, and half intersite distance at — 20 nm, chosen to char- 
acterize J-doped dopants in the insulator with a typical small 
radius and a large intersite separation. Lastly, unless other- 
wise noted, we consider an external magnetic field S = 0.7 
T applied along the device plane, perpendicular to the inter- 
dot axis, providing a Zeeman triplet splitting of 17.5/ieV. This 
field has been used in recent experiments that utilized nuclear 
state preparation to generate Overhauser field gradients,^ii 
and in experiments that studied control sequences to enhance 
qubit coherenceJ^i^ 

In Figs. Q2a) and Q~fb) we examine the couplings depen- 
dence on the polar and azimuthal orientations of the qubit- 
TLF axis, respectively. We fix the TLF polar angle 9t to either 
(vertical traps) or 7r/2 (lateral traps) and look at the coupling 
dependence on 9, <f> for the two 6*r orientations, averaging out 
the remaining two angles. The blue lines in Figs. Q~fa) and 
[Hb) clearly replicate the angular dependence of v\ , given in 
Eq. (O, with minimal dependence on 6*t from higher order 
terms. For v 22 , averaging over the remaining two angles re- 
sults in 

v 22 ~ sin 29 cos <f), 

for the case of 9t — 0, in close agreement with the solid red 
lines in Figs. [Ha), andQlb). For 9t = 7r/2 we have 









- -^j sin 2 9 







vf{(j)) ~ |1 - 3 (cos 20 + sin 20) I . 

Notice that at angles where v 22 changes sign, contributions 
from higher order terms in the multipole expansion become 
important. In addition, the relative strength of these higher 
order terms increases with reduced qubit-TLF distance R, re- 
sulting in a more complicated angular dependence of the cou- 
plings. Both vp and v 1 depend only mildly on the TLF az- 
imuthal angle <pT (not shown). 

Fig- He) shows the R dependence of the couplings at zero 
detuning. Solid lines depict angular-averaged results, whereas 
dashed lines show couplings calculated by fixing the polar 
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Normalized Bias B (T) 



FIG. 1: (color online) qubit-TLF couplings vs. (a) polar angle 8; (b) 
azimuthal angle (f>\ (c) qubit-TLF distance R; (d) TLF Bohr radius 
-Dt; (e) interdot bias normalized to QD confinement and measured 
from the singlet anticrossing point; (f) external magnetic field. In 
all plots blue (red) lines correspond to vp (« 7 ). Left (right) y axis 
in panels (a) and (b) correspond to Vp (u T ). Green lines in panels 
(e) and (f) depict the exchange interaction. Unless otherwise noted 
we consider the singlet anticrossing point (zero detuning), and qubit- 
TLF distance R = 300 nm. Other system parameters are given in 
the main text. 



angle of the TLF-qubit axis to 9 = arccos(L z /i?), where 
L z = 90 nm is the depth of the 2DEG plane below the surface. 
We see that vp ~ R~ 2 for R > 50 nm, and v 7 ~ R~ 3 for 
R > 30 nm, corresponding to the leading contributions in the 
multipole expansion. When 9 is fixed by R, the leading terms 
in both vp and v 1 change sign as R approaches L z (9 — s- 0) 
and higher order terms become dominant. In the rest of the 
paper, having no a-priori knowledge of the TLF location and 
relative orientation, we calculate the couplings by performing 
averaging over all angles^ 

In Fig. Q2d) we show the couplings dependence on the TLF 
site radius D t , where the intersite half separation is varied ac- 
cordingly as at = 4£> t . w 7 shows a linear dependence in D t , 
matching the linear dependence of v 22 on at- On the other 
hand, vp depends very mildly on D t as long as D t < 40 
nm, since its lowest order term that depends on D t and at 
is dipole-quadrupole, which contributes only slightly to vp 
within the range of relevant trap sizes. Fig. [TJe) shows v 7 ,vp, 
and J dependence on the interdot bias point, e, normalized 
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to QD confinement, where detuning is measured from the 
5(1, 1) — 5(0, 2) anticrossing point. The bias dependence is 
qualitatively approximated by the analytical formulas for the 
leading contributions, Eqs. |@), and (0, with 40 



J(s) 



2T 2 



v/e 2 + 4T 2 - e 



Finally, the magnetic field dependence of w 7 , vp, and J, is 
shown in Fig. |TJf). A decrease of more than four orders of 
magnitude in the couplings is observed between B = 100 
mT to B = 1.5 T. We note that the magnetic compression 
of the orbitals reduces the wave function overlap by a factor 
of m 10 within this field range, thereby it alone cannot ac- 
count for these results. In fact, the main reason for the high 
sensitivity of the couplings to the magnetic field is that vary- 
ing B also shifts the bias position of the singlet anticrossing. 
In Fig. [Tff) we have fixed the bias at the anticrossing point 
found for B = 0.7 T. As B increases, the anticrossing shifts 
to a more positive bias, and as a result our fixed bias becomes 
negative, thereby reducing the couplings. In experiments the 
bias position is likely calibrated when B is varied and the high 
sensitivity of the exchange and qubit-TLF couplings would be 
greatly reduced. 

We conclude that the qubit-TLF coupling strength may vary 
by several orders of magnitude and is particularly sensitive to 
the qubit working position (interdot bias) and the TLF dis- 
tance. For concreteness, we consider in the rest of this pa- 
per only 7-type TLFs, given their relative abundance, and will 
henceforth drop the 7 subscript in u 7 to avoid clutter. 



III. TRANSFER MATRIX METHOD 

In this section we outline the formulation for qubit deco- 
herence due to a single classical TLF. 24 We include asymmet- 
ric TLF switching, relevant in cases where the temperature is 
lower or comparable to the TLF level spacing. 

We first write the qubit-TLF Hamiltonian, Eq. (fJJ, in the 
basis of the qubit eigenstates: 



% = Aa' z + v£(t)(smipa' x — cos <pa' z ), 



(6) 



where A = + J 2 , ip = arctan(<$/i/J), and a' x z are 

the Pauli matrices in the rotated frame. In this frame, the qubit 
evolves under a static field in the z axis, with noise in both x 
and z axes. Since we are dealing with a bistable fluctuator, the 
Bloch vector representing the qubit state in the rotated frame 
precesses around either of the two effective fields 



B± = (±v sin ip, 0, A qp v cos ip). 



(7) 



We denote the two TLF states corresponding to B + and B 
as up and down, respectively, and 7+ (t_) is the switching 
rate from state up to down (down to up). Unlike the case 
of pure dephasing (ip = 0), the Bloch vector does not pre- 
cess around the z axis alone but can reach any point on the 
Bloch sphere, resulting in dissipative dynamics. Denoting by 
p(r,t) the probability to reach point r = (x, y, z) on the Bloch 



sphere at time t, we introduce p + (r,t) and p_(r,i) as the 
probabilities to reach the point r when the TLF is in state up 
and down, respectively. 

The equations for p± (r, t) are found to be 

P±(r,t + r) = (1 -l±T)p ± {U ± 1 v,t) + 

7 T rp T (C/ ± 1 r,t), (8) 

where U± = e rB± R rotate the qubit around the B± axes. 
We consider an infinitesimal time step dt, such that j± dt are 
the switching probabilities for the two TLF states within time 
dt. Using the fields B± in Eq. (0, and the rotation matrices 
R = (R x , R y , Rz), 24 we find the master equations: 

p± = [(A =F v cos (p)(yd x -xd y ) ± v smp(zd y -yd z )]p± 
+ (1tPt -1±P±)- (9) 

The expectation values of the Bloch vector components in the 
rotated frame evolving under the fields B± are defined as: 



J drp ± (r,t)r, 



(10) 



and their dynamical equations are found from Eqs. (|9) to be 

x± = {1t x t -1± x ±) - (A T v cos p) y± 

V± = (7t2/t _ 1±V±) + ( A T v cos l f) x ±Tv sin (pz± 

z± = {1t z t - l± z ±) ± vsimpy ± . (11) 

Finally we transform these coupled equations to a set of equa- 
tions for the total expectation values, x — x+ + X- , y = y+ + 



y~,z — z + + z-, and their differences, 5x 

y+ - y-,Sz = 2+ - z_: 



x-,6y 



x = —Ay + v cos tpSy 

y = Ax — v cos ipSx — v sin ipSz 

Sz = — 5jz — 2j5z + v sin tpy 

5x = — Sjx — 2j5x — ASy + v cos tpy 

5y = — 5^/y — 2^/Sy + A5x — vcostpx — vsinipz 

z = vsmipdy, 



(12) 



where we have defined the average switching rate, 7 = (7+ + 
7_)/2, and switching rate difference, £7 = 7+ — 7_. It is 
convenient to write Eqs. ( TT2l in a matrix form 

k = Mik, 

where entries 1,2,6 of the vector 



/ X \ 

y 

Sz 
Sx 
Sy 
V z J 



(13) 
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define the Bloch vector components, and the matrix Mi reads 
/ 




A 


— £7 

■V COS ip 





A 

— L\ 


n 
u 


n 
u 


V COS 







—v sin 


— V COS if 








v sin (/? 


-2 7 








— ($7 


w cos ip 





-27 


-A 





—Sj 





A 


-27 


—ii sin y 











i>sin</5 


/ 
(14) 



We now consider the dynamics of the Bloch vector sub- 
ject to sequences of control tt pulses along an arbitrary axis. 
The application of a 7r-pulse flips the components of the effec- 



tive fields in Eqs. perpendicular to the pulse axis. While 
our dynamical equations are written in the rotated frame, we 
consider 7r-pulses along the original frame axes since these 
are more readily available in current experiments in S — To 
qubits. 7 The resulting effective fields after the application of 
the control pulses read: 

B± x = (A sin 2ip q= v sin <p, 0, — A cos 2<p ± v cos tp) 

= — (±v sin ip, 0, A qp v cos tp) 
B J = (— A sin 2tp ± v sin tp, 0, A cos 2<p =F v cos <p) (15) 



These fields lead to qubit state evolution governed by matrices M|, where j — x,y,z denotes the control pulse axis (in the 
original frame). Mf is found as 



( 





A cos 2(p 








— v cos tp 


° \ 




-A cos 2ip 





vsimp 


v cos ip 





—A sin 2ip 







— v sin ip 


-27 





A sin 2ip 


—(57 




—<57 


— V COS ip 





-2 7 


A cos 







UCOS(£ 


—(57 


—A sin 2ip 


—A cos 2ip 


-27 


vsimp 





A sin 2tp 








—v sin 


/ 



(16) 



M| is found by substituting A — > — A, w 



in Mi, and Mf is found by substituting A 
I 



-A, v 



-v in M 9 2 



The evolution of the qubit state after a sequence of N Ttj- 
pulses can be written as 



k(i)=Z4k(0) 



(17) 



with 



n 



= A«T2i„MlT2i 



2V/2 



TV e odd 



(18) 



+ 1 ]Je M 2 r2 'e Mir2 >-\ AT e even 



where r, is the time interval between the (i — l)th and ith 
pulses. 

When ip = 7r/2, and for symmetric switchings (Sj = 
0), the matrices Mi and Mf become block diagonal and 
the equations for (x,y,8z) decouple from the equations for 
(6x,8y, z)w^ Previously found results for pure dephasing 
(<p = 0) can be obtained more easily by solving a subset of 
Eqs. ( TT2l , as shown in Appendix A. 

For all the results presented below we take the initial state 
of the qubit to lie along the rotated x axis, tilted by an an- 
gle tp in the x-z plane. Such initialization can be realized in 
two steps. First, the initial singlet state is transformed into the 
ground state of the nuclear field, ( 1 5) + 1 To ) ) / y/2, by lowering 
the exchange adiabatically with respect to the nuclear mixing 
time. 4 Second, a magnetic field gradient is ramped slowly so 
that the qubit state follows it, reaching the Hamiltonian eigen- 
state. We remark that this last step is challenging with cur- 
rent experimental techniques that generate nuclear field gradi- 
ent, since the employed dynamical nuclear polarization cycles 



are slower than the electron spin dynamics. Our choice of 
initial state is motivated by the clean context it provides for 
the analysis of the qubit dynamics. In Appendix B, we pro- 
vide formulas for qubit decay, when its initial state is set as 
(\S) + |To))/v2, and discuss its relation to the results given 
in the main text. 

The probabilities of finding the TLF state up or down are 
7 _/(7 + + 7_) and j + / (7+ +7-), respectively, therefore 



r(3 



P±(r,* = 0) = 

7+ +7- 

and the initial qubit's state reads 



<5 (3) (r-x), (19) 



k(0)= 1,0,0 



(5 7 
'2? 



0,0 



(20) 



The ratio between the TLF switching rates is given by the 
Boltzmann factor 



7- = e -AE t /k B T 

7+ 



(21) 



where AE t is the TLF level splitting, taking into account in- 
tersite tunneling. 30 For the TLF parameter range considered in 
this work, AE t < 5/ieV, thus for sample temperature of 100 
mK, we always have <5 7 = 2 7 tanh(A_E t /2fcsT) < 7. 

We now address the choice of rotation axis for the control 
pulses. First, since charge noise is induced along the original z 
axis, we expect that its decohering effects will not be corrected 
by 7r z pulses, which have been extensively employed to reduce 
nuclear induced noise. This is demonstrated in Fig.|2ja) where 
we plot qubit decay under both free induction (FID), and a sin- 
gle 7r 2 pulse (SE). 41 Here we take TLF switching time of 0.1 
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FIG. 2: (color online) Qubit decay under a SE control pulse with 
different rotation axes, (a) ir z pulse (symbols) and FID (solid lines) 
decays for TLF distances of R — 200 nm and R — 300 nm. The 
qubit is at large negative bias, J = 0.02/^eV, and Sh — 0.125/ieV 
(ip < 7r/2). Dashed green lines correspond to the strong coupling 
approximation, Eq. d22t : (b) qubit decay under ir x and n y pulses for 
same parameters as in panel (a); (c) qubit decay under ir x and n y at 
zero detuning, where J = 4.15/xeV 3> Sh — 0.125peV (<p ~ 0), 
for R = 300 nm (v = 1.4/ieV); (d) Same as panel (c) with Sh = 
0.5/j.eV. In all plots the TLF average switching time is 0.1 ms. 

ms and consider a large negative bias such that J = 0.02/ieV. 
For a qubit-TLF distance, R = 200 nm (v = 0.2 neV), we 
have v 3> 7, where the FID decay is dominated by v: 

Xfid ~ e -7 ' (cosvt + — sinutj . (22) 

At R = 300 nm, v = 62peV > 7, and the strong cou- 
pling approximation is compromised [see dashed green lines 
in Fig. [Ha)]. In both cases the tt z pulse does not extend the 
qubit coherence as compared with FID, resulting in dephasing 
times of 22/is, and 95/is, respectively. These relatively short 
dephasing times were likely circumvented in experiments by 
selecting relatively quiet samples. 

The effectiveness of ir x vs. ir y control pulses in mitigat- 
ing charge noise depends on the TLF parameters, as well as 
on the qubit working position. At pure dephasing (Sh = 0), 
both pulse types eliminate noise from TLFs that are static dur- 
ing the time interval between pulses, r, since the qubit state 
evolution depends only on the a z operator [see Eq. (|6]l, and 
Ref.HJfor more details]. At any other working position, n y 
pulses are always superior to ir x pulses, where the relative 
improvement depends on the TLF parameters. The most pro- 
nounced improvement is expected at large negative bias where 
ip < 7r/2, as verified in Fig.|2b). In Fig.|2jc) we tune the in- 



terdot bias to the singlet anticrossing (zero detuning), result- 
ing in J = 4.15/ieV ^> Sh. Since tp « the performance 
of tt x and TT y are similar. However if we stay at the anticross- 
ing and take Sh = 0.5/ieV, increased dephasing is exhibited 
by the ir x pulse [see red dashed line in Fig.|2d)]. Two-qubit 
gates that are performed near zero detuning (to increase cou- 
pling strength) will thus benefit greatly from ir y pulses that 
will allow to increase Sh, shortening gate times. 

Finally we note that it is possible to effectively eliminate 
the exchange interaction so thaty? = tt/2 while retaining sub- 
stantial qubit-TLF coupling. Normally, J and v decrease si- 
multaneously with bias [see Fig. 0I e )L but, as was recently 
shown, at certain interdot bias points, the inter-qubit Coulomb 
couplings balance the internal exchange interaction, resulting 
in a zero effective exchange. 40 For a distance of 400 nm be- 
tween the centers of the two DDs, this zero-exchange posi- 
tion is found at a moderately negative (dimensionless) bias of 
e = —0.05. While this position is typically ideal from the per- 
spective of charge noise, for sufficiently close TLFs, tt x pulses 
will introduce substantial noise. 

The above considerations suggest that n y pulses are more 
effective in mitigating charge noise in several important sce- 
narios. Moreover, we expect tt x pulses to be completely in- 
efficient in correcting nuclear-induced noise, similarly to the 
inadequacy of tt z pulses in correcting charge noise. This will 
play an important role at negative bias, where nuclear noise is 
dominant. While ix v pulses are harder to implement and were 
previously generated with limited fidelity^ it is expected that 
high-fidelity n y pulses will be available in the near future. 42 



IV. ANALYTIC RESULTS FOR PDD AND CPMG PULSE 
SEQUENCES 

In this section we provide exact analytic solutions for the 
qubit dynamics due to its coupling with a single TLF, in the 
limits of weakly (v <C 7) and strongly (v ^> 7) coupled fluc- 
tuators. As we show below, the qubit dynamics is different in 
these two regimes, and the asymptotic behavior in these lim- 
its proves valuable for the interpretation of our results. Our 
analysis follows Ref. |U extending it in three respects: (i) 
we consider an arbitrary qubit working point ip, (ii) we in- 
clude asymmetric TLF switching, and (iii) we find explicitly 
the weights of the different decay rates in the final solution, 
providing further insight into the qubit dynamics, and the en- 
hanced performance of CPMG over PDD control pulses. 

In the following we consider an A-pulse PDD protocol with 
a constant time interval between pulses, r = t/(N+l), where 
t is the total time. (CPMG protocol will be considered in sec- 
tion [IVC]) In light of the discussion ending the previous sec- 
tion, we consider only tt v pulses. In order to calculate the 
qubit dynamics, we first note that after a ir y pulse is applied, 
the matrix that governs the qubit evolution, MM, can be writ- 
ten in terms of M\ in Eq. ( fT4l as: 

Mf = LMiL, (23) 

where L — diag(l, —1,1,1,-1,1). Thus the qubit evolution 
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under the full control sequence is found by 



k(t) = T N+1 k(0), 



(24) 



where (N is odd for PDD) 



(25) 



The decay rates of the Bloch vector components are found in 
terms of the six eigenvalues, xu of the evolution operator T: 



ln|Xi 



(26) 



and the general solution for the decay of the qubit state rotated 
components is 



/(*) = £' 



j' -r 4 t 



(27) 



where is the weight of the ith eigenvalue in the solution 
of the jth component. It should be noted that the decay rates 
defined in Eq. (|26T l and calculated below, are time dependent, 
unlike the more commonly encountered decay coefficients ob- 
tained in Bloch-Redfield theory. Notice also that while we are 
only interested in the dynamics of the three rotated Bloch vec- 
tor components, in the strong coupling case discussed below 
there are nonzero weights for all six eigenvalues of matrix T. 



A. PDD: Weak Fluctuator 

In order to find the eigenvalues Xi of T we need to expo- 
nentiate Mi in Eq. ( TBI . In the case of weak coupling we per- 
form second order (degenerate) perturbation in v/j. Note that 
while Mi is non hermitian due to the unequal switching rates, 
<5 7 ^ 0, we still have Mi = Aff + iM(, where Aff, M[ are 
both hermitian. In this case, the eigenstates of Mi and M\ 
form a bi-orthogonal set and a perturbation theory is readily 
available. 43 The three relevant decay rates are found as: 



H" =( l-|£)3q^5 [(l-A-flOsmV + CcoBV] 



2 .3 i 4 7 2 ; a 2 + 4 7 2 



[(2-Bi-B 2 )sm 2 ip- F 



T-JF 2 + D 2 sin 2 2<p 



(28) 



where the different functions of r are given by: 
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(29) 



Here we used sincf = sinf/f, where f = Ar is the normal- 
ized time interval between pulses. We note that the effect of 
asymmetric TLF switching in this regime amounts to reduc- 
ing all rates by a common factor that becomes appreciable for 
TLF level splitting above 1/ieV '. The three additional eigen- 
values, associated with the dynamics of (5x,Sy,6z), induce 
much faster decay rates, 2 7 — 0(v 2 / 7 2 ), but their weights 
vanish for all three Bloch vector components. Eqs. (128V (|29l 
recover previous results obtained for tp = and tp = ir/2 with 
<5 7 = 0. 24 Specifically, for the case of pure dephasing (ip = 0), 
we have a single decay rate (see also Appendix A): 



r^ = o) = r^ 



i 



tanh 7 r 
7 r 



(30) 



For J = (ip = 7r/2), the equations for x and y decouple 
from z, the latter having no dynamics when the qubit state is 
initially along the x axis, and we obtain two decay rates: 

rr(y = 7r/2)= J V \ J l-A-Bi) 



rs'fc = 71-/2) 



A 2 + 4 7 2 

9 

7 W 

A 2 + 4 7 : 



■(l+A-Ba). 



(31) 



In order to explicitly find the time dependence of the Bloch 
vector components, we diagonalize the matrix T in Eq. d24l >. 
T = SWS^, where W is diagonal matrix, similar to T, whose 
ith element is e~ FiT , and the columns of S are the eigenvec- 
tors of T. Up to second order in v / 7 we can use the unper- 
turbed eigenvectors, and by applying the resulting operator to 
the initial vector, Eq. (f20b . we find the following weights 



J 2,3 



sin — 

2 



= - cos z - 1 ± 



F 



^F 2 + D 2 sin 2 2ip 



■ smr 



u 2,3 = 



1 / F 

—. smr 1 ± 
4 V y/F 2 + D 2 sm 2 2ip / 



w\=Q 



^2,3 
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±— cos — - 



D sin 2tp 
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(32) 
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FIG. 3: (color online) Panels (a)-(c): Weakly-coupled fast TLF 
(I/7 = O.lfis, v = 2.5 neV), A ~ 7 2> v. (a) Decay rates 
vs. time interval between pulses, r, from Eqs. ((28} and j29\ \ (b), 

(c) Bloch vector components under 11 -pulse PDD sequence. Panels 

(d) -(f): Weakly-coupled slow TLF (I/7 = 0.1 ms, v = 8.4 peV), 
A > 7 > t. (d) Decay rates from Eqs. d28t and d29t ; (e) x compo- 
nent under 1 1 -pulse PDD sequence (y and z components are almost 
unaffected by the TLF); (f) Same as panel (e) with J = 0. In pan- 
els (b), (c), (e), and (f), solid blue lines depict Eqs. d27b and d32t . 
and dashed red lines show transfer matrix calculation. In all plots 
J — Sh = 0.125/ieV (ip — 7r/4), except for panel (f) where J = 0. 
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FIG. 4: (color online) Strongly coupled TLF with v = 2.5 neV {R = 
300 nm) and I/7 = 0.1 ms, for Sh = 0.125AteV, satisfying A > 
« 3> 7. (a), (b) Decay rates vs. time interval between pulses, r, 
from Eqs. 03} for J = 5/i = tt/4), and J = (fp = tt/2), 
respectively; (c), (d) qubit decay under 11-pulse PDD sequence for 
J — Sh and J = 0, respectively; (e) Bloch vector y component 
dynamics under 1 1 -pulse PDD sequence (solid blue line) and J = 
(dashed red line); (f) Bloch vector z component dynamics under 1 1- 
pulse PDD sequence for J = Sh (no z dynamics for J = 0). The 
dotted green lines in panels (d) and (e) correspond to Eqs. ((36} and 
{37}. The insets in panels (e) and (f) depict the long time dynamics 
of the transverse Bloch vector components. 



Notice that J2i w f ~ 1> an d Yli w l = Y^i w i = 0- 

In Fig. [3] we illustrate the above results for a working po- 
sition ip — 7r/4 (J = Sh = 0.125/ieV), obtained at nega- 
tive bias e = —0.054. In Fig. Oa) we consider a fast TLF 
with average switching time of O.l/xs, and coupling strength 
v = 2.5 neV (R = 300 nm), such that A ~ 7 > v. All 
three decay rates in Eqs. (128} are comparable and contribute 
to the decay of all three Bloch vector components shown in 
Figs-Ob) and (c), where we consider the qubit dynamics un- 
der an 11-pulse PDD protocol. In Figs. Od)-(f) we consider 
a slow TLF with an average switching time of 0.1 ms. To 
stay in the weak coupling regime, we take R = 2000 nm, re- 
sulting in v = 8.4 peV, so that A S> 7 S> v. In this case 
the function C in Eqs. d29l ) dominates and we obtain two dis- 
tinct decay rates: Tf ~ r| ~ f (1 - ta " h T 7T ) cos 2 ip, and 

T l 2 fa "^gi- sin 2 ip, the latter being independent of r. Whereas 
Ff, r™ 3> T2, the weight is vanishingly small, and over- 
all the qubit state decays with a single (faster) rate. We also 
observe that in this limit there is almost no dissipative dynam- 
ics since w\ ~ 0, and w\ — —w^ for the y component, and 



W23 ~ (7/ A) 2 <C 1 for the z component. Finally, we note 
that the leading term in the faster decay rate, Tf, is eliminated 
at (p = tt/2, and the qubit decay is governed by the much 
slower rate T™, as well as by higher order corrections in 7/A, 
and v/j to Tf and Tf. This case is shown in Fig.[3jf), where 
we set J = and obtain more than 5 orders-of-magnitude 
increase in coherence time. 



B. PDD: Strong Fluctuator 

In the limit of strong coupling, 7 <C v, it is sufficient to 
perform first order perturbation in j/v. Here the unperturbed 
eigenvalues of T form two degenerate subspaces of dim = 2, 
and dim = 4. The analytic results are complicated and we 
present here only the leading term (third order) in the expan- 
sion in v/A, in the case of symmetric switching (#7 = 0). 
The results below are thus valid at negative detunings, where 
v <C A (for a fixed Sh = 0.125/ieV), but will be applicable to 
an extended bias regime if larger Sh is considered. The three 
decay rates associated predominantly with the Bloch vector 
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FIG. 5: (color online) Strongly coupled TLF with v — 2.5 neV 
(R = 300 nm) and I/7 = 0.1 ms, for Sh = O.Ol^teV, satisfying 
A > v ^> 7. (a), (b) Decay rates vs. time interval between pulses, 
r, from Eqs. &33\ for J = Sh (ip — 7r/4), and J = (<p = 7r/2), re- 
spectively; (c) qubit decay under 1 1-pulse PDD sequence for J — Sh 
and J = 0; (d) y and 2 components dynamics under 1 1-pulse PDD 
sequence (no z dynamics for J = 0). Notice that for comparison 
with Fig. [4] we kept the same coupling strength, although normally 
it would reduce with J. One can tune J without changing the bias 
(thus keeping v fixed) by introducing a second double dot. 



components are found as: 



rf = yAsin z <p + — cos 2 ip 



12,3 = sin 2 <p-F±^F 2 +D 2 sin 2 2^ j. (33 1 

where the different functions are given by 

A = 1 — sincf 



B = 2 1 - sine 



,r 
2 

r f 

D = cos sine— 

2 2 

F = - (a sin 2 93 - ^- cos 2 95 



(34) 



We consider strongly coupled TLF in two scenarios: A ^> 
v ^> 7, and A > v 7, shown in Figs.|4]and|5] respectively. 

Examining the asymptotic behavior of the decay rates at 
f <C 1, we have A = B = —2D = f 2 /6. As a result, at short 
time intervals, the qubit dynamics is governed by a single de- 
cay rate: rf w T| w jv 2 t 2 /6, depicted by the solid blue and 
dotted red lines in Fig. EJa). T| scales as r 4 and is consider- 
ably smaller. At the limit of f ^ 1 we have A, B, D <C f, 



thus rf « r| 



-yv 2 r 2 cos 2 ip/6 > T§ 



27U 2 sin (p /A 2 



Note that for pure dephasing, both asymptotes coincide and 



r% = o) 



7?r 



V 2 /6 is the single solution to Eqs. (1331 . 



We conclude that when u <C A, unless ^ = 7r/2, the qubit 



dynamics is governed at all times by a single decay rate that 
scales as r 2 . In contrast, when J = 0, the cos 2 (f term is 
removed and Tf 2 decrease substantially in the large f limit, 
becoming comparable to Tg. In this case the rates read: 



r s 1 (^-7r/2) = ^(l- sincf) 



1 - 2sinc 2 - 



sincf j 



(35) 



both contribute to the qubit decay, and their magnitude in 
the large f limit is curbed at a value r yv 2 /A 2 , as shown in 
Fig. H|b). As a result, decay time is nearly three orders of 
magnitude longer, as compared with the J = 5h case. This 
is shown in Figs.Uc), an d (d) for an 1 1 pulse PDD sequence. 
We note that in this case, the equations for (x, y, 5z) decou- 
ple from those for (Sx, Sy, z), and T| is associated solely with 
the z component, having no weight in the decay of our initial 
state. 

In order to fully capture the qubit dynamics, we calcu- 
late the weights of the different rates in Eq. (|33l . As in the 
weak coupling case, the three additional eigenvalues, associ- 
ated predominantly with the dynamics of (6x, 6y, Sz), induce 
much faster decay rates: 27 — 0(jv 2 / A 2 ), but in this case 
they have a small yet nonzero weight that contributes to the 
decay of the Bloch vector. We demonstrate this behavior for 
ip = it/ 2, where the analytic formulas are simpler due to the 
block diagonal form of Mi in Eq. ( TBi i. Here we only need to 
consider the three eigenvalues of the block (x, y, 5z), where 
two decay rates are given in Eqs. (l35t . and the third rate is 
T3 27, associated predominantly with 8z. Considering 
the perturbed eigenvectors of T to second order in v/A, the 
weights for the x and y components (z has no dynamics) are 
calculates to be: 




C0S 2 + A^ 



— A 2 
A 2A 



T 



A 



(36) 



and 



,,y — 



2 Sinf+ 4A^ 



A+-(B-2) 



2 2 4A 2 



(B-2)[--A 



A 



;A(B-2), 



(37) 



satisfying J^. wf — 1, and ^\ w\ = 0. The dashed red 
lines in Figs.|4fd) and (e) are plotted from Eq. d2~7t . using the 
weights in Eqs. P61) . and (I371 l. whereas the dynamics of the z 
component, present only for J ^ 0, is shown in Fig. [4{f). 

Turning to the case of strongly coupled TLF with A > v 3> 
7, an analysis similar to the one above qualitatively captures 
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the dynamics, although the analytic decay rates in Eqs. (l33l l 
are accurate only to third order in v/A. To illustrate this 
regime we consider Sh = J = O.Ol^teV, and take J = Sh 
or J = 0. The asymptotic behavior at the short and long r 
limits remains as before, but the crossover from the r 4 short- 
time scaling to the constant long-time behavior of T§ occurs at 
a much larger r, as seen in Fig.[5Ja). Since the qubit dynamics 
are governed by a single decay rate: Tf = T 2 , which is inde- 
pendent of A, the resulting qubit decay depicted by the solid 
blue line in Fig. [5jc) is largely unchanged. In contrast, for 
J = 0, the two contributing decay rates scale like 1/A 2 and 
are almost two-orders-of-magnitude larger than the rates plot- 
ted in Fig.Hfb). As a result, a dramatic decrease in coherence 
time is shown by the dashed red line in Fig. 0c), as compared 
with the A 3> v 3> 7 case. Finally, the dissipative dynamics 
in the y and z components are much more pronounced for the 
A > v case, as seen in Fig. [2d). 



C. CPMG 

Here we extend the above analysis to treat the CPMG proto- 
col, where 73 = t/N for 2 < j < N and t x = t n+ i = t/2N. 
The evolution of the qubit state under TV-pulse CPMG se- 
quence can be written as: 



k(t) 



T 1/2 T N - l T 1/2 ,k(0), N e odd 
LT 1/2 T N - 1 T 1/2 k(0), N <= even, 



(38) 



where Ti/ 2 is defined by Eq. d25l l with r — > r/2. Focusing 
on the x component of the Bloch vector, we consider first a 
weakly coupled TLF (v <C 7), in the case of J = 0. Here 
the z component is again decoupled from x and y, and we can 
consider only the 3x3 block in T corresponding to (x, y, Sz). 
Furthermore, since there is no weight for the fast ~ 27 decay 
rate, predominantly associated with Sz, in the weak coupling 
case, we can work within the 2x2 subspace of (x, y). 

The similarity transformation that diagonalizes our subset 
of T is: 

and we can write Eqs. (|38l for the (x, y) components as: 





x(t) 

y(t) 



LS 



Xi u \ of 

jV-1 1 D 
A2 







(39) 



where Xi are the two eigenvalues of T that are associated 
with the decay rates in Eqs. (fjlj [recall Eq. (|26|il, and L = 
diag(l,±l), with the upper (lower) sign corresponding to 
odd (even) number of pulses. Note that x(t) is unaffected by 
L, and we obtain an 'even-odd effect' only for the y com- 
ponent. We find that x(t) decays with a single rate, TJf, 
given in Eqs. OTb . The well known superior performance 
of the CPMG protocol over PDD can thus be explained as 
follows. For 7 <C A, Eqs. ( BTT l can be approximated by 
w ^(l-sincf), andT™ « ^(l + sincf-2sinc 2 f/2). 



An efficient noise suppression is thus achieved when the time 
interval between pulses satisfies f <C 1. At these short time 
intervals, Tf ~ 7« 2 t 2 > - 7u 2 A 2 r 4 . Since only T 2 is 
present in the CPMG pulse sequence, it is more effective than 
PDD at the 99 = 7r/2 point. 

Similar calculation can be carried for a general working po- 
sition. In the general case, we need to include the z compo- 
nent that is now coupled to x and y. We find again that the x 
component decays with a single rate given by: 



"fV 



A 2 + 4 7 2 



[(1 + A- B 2 ) sin 2 (p + C cos 2 if] , (40) 



where A, B 2 , and C were given in Eqs. (|29l . The cos 2 ip term 
in Eq. (l40b dominates at short r, and the advantage of CPMG 
over PDD is largely eliminated at ip 7^ tt/2. 

Turning to the case of a strongly coupled TLF {v ^> 7), 
we include the small weight of the fast decaying (5x, 5y, Sz) 
components. Considering only the ip = tt/2 case, we per- 
form a calculation similar to the above within the subspace 
(x, y, Sz), with Tf 2 given in Eqs. ( f35l > and the additional 
small weighted rate of 27. As in the weak coupling case, we 
find that the weight of Tf is vanished and the resulting qubit 
decay reads: 



2 \ e -nt 



+ ^A 2 e- 2 ^\ (41) 



where A was given in Eqs. d34l l. While the weight of the sec- 
ond term is small, its fast decay rate induces a sizable effect 
on the qubit decay. 



V. TLF-INDUCED DECOHERENCE UNDER GENERAL 
PULSE SEQUENCE: RESULTS 

In this section we present qubit coherence times due to its 
coupling to a single TLF. With limited knowledge of the TLF 
characteristics, we wish to map a wide range of TLF switch- 
ing rates and coupling strengths, at various qubit bias points. 
For brevity we focus only on the qubit decay time, T 2 , which 
we define as the time it takes the x component of the Bloch 
vector of a qubit initially prepared along the x axis to reach 
50% of its initial value. We remark that it is often desirable 
to study the initial qubit dynamics, which may be different for 
two pulse protocols that have similar T 2 . In addition, T 2 does 
not capture the dissipative dynamics that leads to nonzero y 
and z components (see Figs. [3}|5j. Nevertheless, the results 
presented below give good indication for the qubit dynamics 
in the various regimes. 

As a benchmark case, we consider a 10-pulse CPMG se- 
quence. Fig. [6] shows T 2 times vs. TLF average switching 
times for R = 300 nm, at three bias points. We identify three 
coupling strength regimes at which very different dynamics 
occur, as explained below: (i) very strong coupling, (ii) strong 
coupling, and (iii) weak coupling, marked by black, green, 
and red dashed lines, respectively. Note that by weak and 
strong we refer only to the ratio v/j. In Fig. |6|a) we con- 
sider zero detuning, where J = 4.15/LteV 3> Sh (ip ~ 0), 
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lowing dephasing times: 



y 1 (ms) 

FIG. 6: (color online) Qubit decay times vs. TLF average switching 
times under a 10-pulse CPMG sequence, for a qubit-TLF distance 
R = 300 nm. (a) Zero detuning (J = 4.15^eV, v = 1.4/xeV); (b) 
moderate negative detuning (J = 0.125/ieV, v — 2.5 neV); and (c) 
large negative detuning (J = 0.02/ieV, v = 62 peV). In all plots 
<5/i = 0.125/ieV. The black, green, and red dashed lines correspond 
to very strong, strong, and weak coupling regimes, respectively. The 
vertical dotted lines in panels (b) and (c) mark v — 7 (in panel (a) it 
is not within range). 



v = 1.4/ieV, and the qubit is most susceptible to charge noise. 
Here we are in the very strong coupling regime for the entire 
range of TLF switching rates. Since v/A is not small, it is 
more straightforward to consult the results for pure dephas - 
ing given in Appendix A. In this quasi static regime, for our 
case of 10 pulses, the r values relevant for the qubit T% times 
satisfy r ~ I/7 ^> 1/v for almost the entire range. The 
large r limit of Eq. dAl 5b results in T s ± « 7 for both rates, 
thus the qubit approximately decays with a single rate, 7, and 
T 2 =ln2/ 7 . 

Turning to the case of large negative detuning, where J = 
0.02^eV <$:Sh(tp< tt/2) and v = 62 peV, Fig.^c) includes 
all three coupling strength regimes. Focusing first on the weak 
coupling (left side of the vertical dotted line), the qubit decays 
with a single rate given by Eq. d40t , For A ^> 7, satisfied 
here, the cos 2 ip term dominates, as long as we are not too 
close to the ip = tt/2 point [see Eqs. j29l l. and the decay rate 
reads: 



tanIi7T 



cos If. 



(42) 



27 V 7T 
The long- and short-time asymptotes of Eq. (l42l yield the fol 
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rpSt 

1 2 — 



2 7 In 2 

V 2 COS 2 (f ' 



6A 2 In 2 
7« 2 cos 2 ip 



1/3 



7T > 1, 
7T < 1 



(43) 
(44) 



where N is the number of pulses. Observing the solid ma- 
genta line in Fig.[6|c) that depicts Eq. d42l . and the dashed red 
and green lines, corresponding to its asymptotes, we notice 
that Eq. d42b works well into the strong coupling regime. This 
can be explained by taking the large r limit of Eqs. d33l . where 
we have found a single decay rate: V s = jv 2 t 2 cos 2 <p/6, 
which is identical to the large r limit of Eq. d42l . For very 
weak coupling [left end of Fig. [6{c)], the approximate re- 
sult of Eq. (1421 breaks down since A ^> 7 no longer holds. 
We note that although J <C Sh, the dynamics shown here is 
very different from the case of J = 0. In the latter case, T 2 
times increase substantially, due to much smaller rates given 
by Eq. d40b . To realize this limit one needs to work at a very 
large negative bias, or to utilize inter-qubit couplings that can 
turn off the exchange completely. Fig.[6jb) shows decay times 
for moderate negative bias where J = Sh, exhibiting the com- 
plex dynamics discussed above. 

In Fig. [7] we show T2 times vs. qubit-TLF distance, for 
a TLF average switching time of 0.1 ms, at the same three 
bias points. At zero detuning we are again at the very strong 
coupling regime throughout the considered range, resulting in 
a constant T% = In 2/7 = 69.3/is, as shown in Fig. [7Ja). At 
large negative detuning, Fig. [3c) exhibits the weak and strong 
coupling regimes captured by Eqs. ( 1431 , and (l44l . Here, the 
relevant r values scale as 1/v. To the left of the vertical dotted 
line, marking v — 7, we have strong coupling, thus jt <£L 1, 
and Eq. d44l applies. Since the dominant term in the qubit- 
TLF coupling scales as R~ 3 [see Fig. [He) and discussion 
therein], T 2 scales like R 2 in this regime. To the right of the 
v = 7 line we have 77 > 1, and Eq. d43l applies, thus T 2 
scales as R 6 . The T 2 times in the intermediate bias position 
shown in Fig.[7fb), are understood by the same arguments all 
the way down to R as 250 nm, below which, the very strong 
coupling, quasi-static regime applies, and T 2 is roughly con- 
stant. The increase in T 2 times below R = 100 nm is not 
explained within our analytic results, which do not apply to 
this v 3> A regime. 

The discussion above suggests that prolonging the qubit co- 
herence time by reducing the time interval between pulses 
(or equivalently increasing the number of pulses) is effective 
only when Eq. ( l44b holds. Within this regime, marked by the 
overlap of the dashed green lines with the full calculation in 
Figs-121 and[7] decay time scales with the number of pulses as 
N 2 / 3 . This power law was predicted for pure dephasing due 
to random telegraph noise in Ref. There it was explained 
by observing that as N increases, the non-Gaussian noise at- 
tributes are suppressed by the DD pulses. The A 2 / 3 scaling 
was then found by considering Gaussian noise with spectral 
density having a soft cutoff (e.g. w~ 2 ). 23 Furthermore, the 
same power law was observed in experiments on spin qubits in 
nitrogen vacancy centers in diamond, where both CPMG and 
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FIG. 7: (color online) Qubit decay times vs. qubit-TLF distance 
under a 10-pulse CPMG sequence, for a TLF average switching time 
of 0.1 ms. (a) Zero detuning (J = 4.15^ieV); (b) moderate negative 
detuning (J = 0.125/ieV); and (c) large negative detuning (J = 
0.02/ieV). In all plots Sh — 0.125/ieV. The black, green, and red 
dashed lines correspond to very strong, strong, and weak coupling 
regimes, respectively. The vertical dotted lines in panels (b) and (c) 
mark v = 7. 



double-axis (XY) DD pulse sequences were applied. 44 In this 
system, the dephasing of the central spin is induced by its cou- 
pling to a bath of spins with dipolar intra-bath coupling. By 
identifying the noise as a Gaussian and Markovian Ornstein- 
Uhlenbeck process, the authors were able to explain the TV 2 / 3 
scaling using the arguments given in Ref. |23j 

We find estimates for the number of pulses required to 
crossover into the large-iV regime from either side. At the 
very strong coupling regime (dashed black lines), T2 — 
In 2/7, therefore to avoid the large r limit of Eq. ( |A15t , we 
require vt w (v/7) In 2/N < 1, or N > (u/7) In 2. At the 
weak to strong coupling crossover we have T2 = In2/T st , 
given by Eq. d44i i. and to satisfy the condition 7T < 1, we 
need N > 6 In 2(7/1; cos ip) 2 . To summarize, the minimum 
number of pulses needed to enter the regime at which T2 
scales as N 2 / 3 is estimated by: 



N > 



v\n2 61n2 Y 



7 



cos 2 (p v 2 



(45) 



These considerations are demonstrated in Fig. [8] where we 
plot T2 times vs. the number of CPMG pulses for a qubit- 
TLF distance of R = 300 nm, at the above three bias points. 
The dashed lines correspond to the power law: T2 oc N 2 / 3 . 
The vertical dotted lines mark the minimum number of pulses 
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FIG. 8: (color online) Qubit decay times vs. number of CPMG pulses 
for R — 300 nm. (a) Zero detuning (J = 4.15/ieV); (b) moderate 
negative detuning (J = 0.125/^eV); and (c) large negative detuning 
(J = 0.02/ieV). In all plots Sh = 0.125fj>eV. The dashed lines cor- 
respond to Eq. d44t . The vertical dotted black (red) lines mark the 
crossover from very strong (weak) to strong coupling regime. 



necessary to enter the above regime, where black and red lines 
correspond to the left and right conditions in Eq. (l45t . respec- 
tively. While the N 2 / 3 power law is in effect only for N satis- 
fying Eq. ((45), considerable increase in coherence time is still 
obtained for smaller number of pulses, to the left of the ver- 
tical dotted lines. It is interesting to note that the 2/3 power 
law is close to the value 0.72, recently found for even num- 
ber of CPMG pulses in a setting that was likely dominated 
by nuclear-induced dephasing. 13 Finally, we comment that we 
have assumed that our control pulses are ideal, in that they 
are zero-width tt pulses applied exactly along the y axis. In 
reality, control pulses have errors in both their rotation angle 
and axis, which may introduce more noise than they can ef- 
fectively remove^i 4 ^ This is particularly relevant for large- N 
sequences, thus the estimates in Fig.[8]should not be taken too 
seriously above N f=s 50. 

Next we compare T2 times under several widely used DD 
protocols. The concatenated DD sequence (CDD) is defined 
recursively, 46 with the Zth order of concatenation given by: 

CCQ(i) = CDD i _ 1 (i/2) - tt- CDD;_i(t/2), (46) 

where CDDo is free induction, and the middle n pulse is in- 
cluded only for odd /. The number of pulses in the Zth level 
concatenation is found as: 



N, 




l) , l e odd 
2J , I € even. 



(47) 
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FIG. 9: (color online) Qubit decay times vs. number of pulses for 
R — 300 nm under various protocols, (a) Zero detuning (J = 
4.15/ieV), TLF average switching time of 1/is; (b) moderate neg- 
ative detuning (J = 0.125/ieV), TLF average switching time of 0.1 
ms; and (c) large negative detuning (J = 0.02/xeV), TLF average 
switching time of 1 ms. Panels (d), (e), and (f) show decay times to 
95% of initial value, for same parameters as in panels (a), (b), and 
(c), respectively. In all plots Sh — 0.125/ieV. CDD2 and 2-pulse 
UDD are the same as 2-pulse CPMG, and PDD has only odd number 
of pulses. 



The Uhrig DD sequence (UDD) is defined by pulse times: 



U = t sin 



2(JV + 1) 



1< i < N 



(48) 



where the time interval between the (i — l)th and ith pulses 
is n = U — ti-i (t = 0, £at+i = t). A1 We have not at- 
tempted to obtain analytical results for the decay rates of pulse 
sequences with unequal time intervals between pulses. In 
Fig. |9]we compare the performance of PDD, CPMG, CDD, 
and UDD sequences at the above three bias points. At zero 
detuning, Fig.|9ja) shows no advantage in increasing the num- 
ber of control pulses up to TV = 50, for any of the proto- 
cols, even though we are considering a relatively fast TLF 



(7" 



Ins). At this strong coupling limit, the qubit state 



decay is characterized by plateaus (see, e.g., Fig. [Tol l that re- 
sult in the observed oscillations in the dependence of T2 on 
N. A better performance measure may be decay time to 95% 
of initial value, plotted in Fig. |9|d), where coherence times 



are consistently improved with the number of pulses. Decay 
times at moderate- and large-negative bias points are shown in 
Figs.|9lb) and (c), respectively, with their respective 5% drop 
times given in Figs. [9je) and (f). Throughout the wide TLF 
parameter range considered, and at all qubit bias points, the 
CPMG protocol exhibits superior performance in mitigating 
charge noise. These results agree and complement previous 
observations of the superiority of CPMG over UDD in fighting 
noise with power-law high-frequency taili2ii£ We also note 
that all pulse protocols roughly follow the N 2 ^ 3 power law 
over a wide range of TLF parameters and qubit working posi- 
tions. 



VI. CONCLUSION 

In this work we have studied the dynamics of a singlet- 
triplet spin qubit afflicted by a charge fluctuator, under dy- 
namical decoupling control pulses. We have presented a the- 
ory that predicts rich dynamics governed by the qubit-TLF 
coupling strength and the TLF switching rates, as well as by 
the qubit working position. For a relatively small fixed mag- 
netic field gradient, the exchange interaction at zero detuning 
results in pure dephasing that is predominantly dependant on 
the TLF parameters. In contrast, at negative bias points the 
qubit dynamics becomes dissipative and the effectiveness of 
the DD pulses depends on both TLF and qubit characteristics. 

We have demonstrated that tt v pulses are preferable over 
ti x pulses in several scenarios, and moreover we expect that 
TT y pulses will be superior in eliminating nuclear fluctuations. 
Finding analytical formulas for the qubit decay rates in the 
limits of weak and strong coupling for PDD and CPMG se- 
quences, enabled us to explain our results for qubit coherence 
times as functions of TLF distance and switching rate at var- 
ious bias points. Over a large range of system parameters, 
coherence times follow a power law T% oc N 2 / 3 , where N is 
the number of pulses. In addition, comparing the performance 
of several pulse protocols, we found that CPMG is the most 
effective protocol to eliminate charge noise. While we have 
presented specific results for two-spin qubits in gate-defined 
GaAs QDs, the formulation presented in this paper should be 
relevant to a wide variety of systems afflicted by charge fluc- 
tuators. A natural extension of our work will include pulse- 
error analysis^^ In particular, the tolerance of various pulse 
sequences to accumulation of systematic and random pulse 
errors should be evaluated in the context of charge noise. 

As current experiments move to more complicated QD 
structures, such as two coupled double QDs, 1011 and three- 
spin qubits,^£ we expect that charge noise will play an in- 
creasing role in the system decoherence. There are two main 
aspects needed to be considered when analyzing charge noise 
in larger devices. First, as the sample size increases, larger 
number of active charge traps can couple to the qubit(s). It 
is well known that a collection of weakly-coupled TLFs with 
a broad distribution of switching rates can lead to 1/ f noise 
spectrum. On the other hand, qubit decoherence due to small 
mesoscopic ensembles can be dominated by the strongest fluc- 
tuator(s), and the generated noise is non-Gaussian with large 
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variability between samples. The system dynamics can thus 
be quite different when crossing over from small to larger de- 
vices. A second point is that coupled qubits may provide ac- 
cess to working points that are unavailable for a single qubit. 
As the interdot bias controls simultaneously the exchange in- 
teraction and TLF coupling strength, the results of section [V] 
were restricted to the available subset of system parameters. 
With two or more qubits, inter-qubit couplings provide an ad- 
ditional handle over J, which is unrelated to the TLF cou- 
pling. In particular, optimal bias points at which J = 0, be- 
come available, where a very different dynamics is predicted. 

Finally, experimental measurement of the qubit coherence 
time under DD sequences can provide insight into the spectral 
characteristics of the dominant noise processes.^^ The dis- 
tinctive dynamics found for weak and strong fluctuators un- 
der DD, can therefore facilitate the characterization of charge 
fluctuators in solid-state devices. Moreover, by comparing our 
calculated decay rates for a general working point, where dis- 
sipative dynamics occur, with the results of a theory based on 
a Gaussian and Markovian noise, one can find both dephasing 
and relaxation functions, and determine the validity range of a 
Gaussian noise theory for qubit-TLF coupling under dynami- 
cal decoupling. 
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APPENDIX A: PURE DEPHASING 

For pure dephasing (93 = 0), Eq. (0 reads B± = (0, 0, J=F 
v). This case is realized when Sh = and has been studied 
extensively. We outline its solution here to provide context to 
the general case studied in the main text. The controlled z ro- 
tation generated by J is eliminated in all cases except for free 
induction, thus we disregard J with the understanding that a 
qubit initially prepared along the x axis acquires only a ran- 
dom phase. Since we need to consider only the x component 
dynamics, it is more straightforward to evaluate the signal de- 
cay by: 



x (t)= / d<t> P (<p,ty*. 



(Al) 



Dividing the probability distribution to p± (c/>, t) to accumulate 
phase <p while the TLF is in the up or down state, we can write 
a set of coupled equations for x±, analogous to Eqs. (|9]): 



X± = ~7±X± + 7tXt T w X±- 



(A2) 



These can be converted into equations for \ = X+ + X— > an d 

Sx = x+ - X-- 



X = -iv6\ 
Sx = -&1X - 2 7^X - ivx, 



which, in turn, are transformed into a second order equation 
for x- 

X + 27X + v 2 x — ivSjx = 0, (A4) 

with initial conditions %(< = 0) = 1, and %(t = 0) = 
ivt&nh(AE t /2k B T). The general solution to Eq. (lA4t 

j 3 2 1.22,24 

X (t) = a+e-^ 1 -^ + a-e-^ 1+ ^\ (A5) 

where 



1- 



1%v , f AE T 
— tanh — — — 
7 \2k B T 



(^ fl ±l)(l±^j) 
2/x 



(A6) 



and /ifl, fij are the real and imaginary parts of fi. The same 
result is obtained from the more general formulation presented 
in section [HI] For pure dephasing with qubit initially along 
the x axis, z± (t) = 0, and the remaining four equations in 
Eqs. (fT2] > decouple into two sets (x, Sy), and (Sx, y), with a 
solution given by appropriately rotating Eqs. (!A5t - (lA6l . 

In order to consider the effects of a sequence of ir pulses 
we first note that in the case of pure dephasing, tt x and n y 
pulses flip the qubit state in the same way, since the qubit 
evolution involves only the a z operator (Sh — 0). Repeated 
applications of tt pulses alternately change the sign of the last 
term in Eq. (IA4b . resulting in a solution similar to Eq. ( |A5) 
with /j, — > n*. A solution for the qubit signal decay after the 
application of a general sequence of tt pulses can be found by 
stitching the solutions between pulses. This is done using a 
transfer matrix approach analogous to Eqs. ([T7l - ([T8V For a 
sequence of N pulses we find: 



and 
where 



4 +1 
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j E odd 



j G even 



-7(fJ+l)TJ / p 2 lt lTj (-1 _ 



(Al) 



(A8) 



(A9) 



Tj is the time interval between pulses j and j + 1, A = A*, 
and a\. are given by Eq. (|A6t . The final solution for the qubit 
signal decay after N control pulses reads 



Xdd 



(t) = a^e"^ 1 -")' + a N+1 e-^ 1 +^ t , 



(A10) 



for even N, and a similar solution with ji — > fi* for odd N. 
As an example we retrieve below the known result for SE, as 
well as the signal decay in the case of two-pulse CPMG 
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Xcp 2 ( t ) = 



-it 



4(/4-l)(M! + l)sinh 



£>n±l)(l±iM/) [(m| - lje^w* 



+ (1 + A*j)c 



2\ ±7/i H t] 



(A12) 



For AT pulse PDD or CPMG sequences, analytical results 
are presented in section [TV] For PDD with symmetric TLF 
switching, it can be shown that the qubit dephasing is given 
by two exponentials, similarly to Eq. dAlOt . These expo- 
nentials are the eigenvalues of the matrix A e (r)A (r), where 
t = tj (N + 1). They are found as;2i 



X± = 



sinh 7/iT ± W cosh 7 /it ^ 



(A13) 



In the limit of weak coupling, v <C 7, the \+ eigenvalue dom- 
inates, and we obtain a single decay rate by using Eq. 



r w = — 1 

27 



tanh 7r 
7r 



(A14) 



This result was found in the main text, by taking ip = in 

Eqs. (EU- and <|29]l. 

In the limit of strong coupling, v 7, Eq. JA13t results in 
two oscillating rates: 



7 (1 ± sinc^T) . 



(A15) 



which contribute to the qubit dephasing with corresponding 
weights: 



w± 



I, 

-(1 ± cosvt) 



(A16) 



The short t limit of Eq. ( lA15b can be retrieved from Eqs. d3Jt 
and d34l i in the main text, where we presented the strong cou- 
pling results for the case of v -^i A. 

Fig. [TO] shows signal decay for several DD sequences, cal- 
culated for a TLF with an average switching time of 1 ms, 
and coupling strength of v = 0.1 neV. The solid blue and 
red lines correspond to SE and 2-pulse CPMG given by 
Eqs. ( IA1 11 1 and ( 1A12| |. respectively. The green line, depict- 
ing 10-pulse CPMG, was generated using the transfer-matrix 
method, Eqs. ( lA7t -( |AT0l i. To verify our analytical results, we 
performed a numerical simulation of the TLF random switch- 
ing, using a Poisson process with time constants j±. The 
symbols in Fig. [TTJ depict an average over 1000 realizations 
of the random switching, with appropriately weighted initial 
TLF states. We have carried similar simulations to verify our 
results for other DD sequences at general working positions, 
presented in the main text. 



APPENDIX B: Effects of frame rotation on qubit dynamics 

The analysis presented in the main text was applied to find 
the dynamics of the qubit state components in the rotated 
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FIG. 10: (color online) Signal decay for ip = under SE (blue line, 
circles), 2-pulse CPMG (red line, triangles), and 10-pulse CPMG 
(green line, squares). Solid lines show analytic results and symbols 
depict averages and standard deviations of 1000 numerical realiza- 
tions of the TLF random switching process. In all plots 5h — 0, 
v = 0.1 neV, average TLF switching time is 7 = 1 ms, and TLF 
level splitting is AEt = 0.45/ieV. 

frame. The initial qubit state defined in Eq. ( f20l > is thus ap- 
propriately set to lie along the rotated x axis. Current experi- 
mental techniques in QD S — Tq spin qubits have limited ac- 
cessibility to this state. The purpose of this appendix is to pro- 
vide a more direct connection to currently investigated S — Tq 
qubits by including the effects of the frame rotation. In this 
appendix only the Bloch components are defined in the origi- 
nal frame, thus the x axis refers to the state (\S) + \T ))/y/2. 
Rotated quantities stated in the main text are denoted here 
with a prime. We note that for pure dephasing (ip = 0) the 
two frames coincide and all our previous results are intact. 
The qubit dynamics at the anticrossing point, where J 5h 
and charge noise dominates, is thus adequately captured in the 
main text. 

Eq. (O, from which our dynamical equations are derived is 
defined in a frame rotated by an angle ip in the x — z plane. 
Vectors and matrices are transformed to this frame by: 

k' = Ok (Bl) 
M' = OMO- 1 , (B2) 



where 



O-- 



/ cos ip 
1 



\ sin ip 






cos ip sin ip 
- sin ip cos ip 
1 




- sin ip \ 





COS(,9 / 



(B3) 



In the following we consider a qubit initially prepared along 
the x axis: \tp(t = 0)) = {\S} + |T >)/%/2. In the rotated 
frame the initial state replacing Eq.(l2"0t reads: 



k'(0) = ( cos <p, 0, — simp 



5~f 5j , 

— ,-cos^— ,0,sin<£ 



(B4) 



To find the qubit's evolution under iV-pulse PDD control se- 
quence we rotate Eq. d24T i back to the original frame: 



k(t) =O- 1 T N+1 Ok(0), 



(B5) 
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where T is defined in Eq. (1251 1. The decay rates, I\, found 
in Sec. IV from the eigenvalues of T, are unaffected by the 
above rotation, and only their weights in the solution for the 
Bloch components are modified. In order to find the trans- 
formed weights we use the eigenvectors of T, as described in 
Sec. IVA. 

For the case of weak coupling the weights of the three decay 
rates, Eqs. d28l l. in the three (non-rotated) Bloch components 
are found as: 



cos ip sin 
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— sin 2os 

4 ^ 



it 



y/F 2 + D 2 sin 2 2^ 



cos z - 1 ± 

2 V y/F 2 + D 2 sin 2 2<p y 



(B6) 



where f , L>, and F are defined in Eqs. (|29] i, and £\ wf = 1, 
and^V w\ — J2i w i = 0- Note that the rotation matrix mixes 
x and 2 components only when ip 7^ 0,7r/2, thus dissipative 
dynamics is minimal when Sh ^> J (large negative detuning) 
or J Sh (at the singlet anticrossing). 

For strong coupling there is a nonzero weight for the fast 
27 decay rate predominantly associated with the dynamics of 



(Sx, Sy, Sz), in addition to the contributions from the three 
rates given in Eqs. ( |33l >. For ip = ir/2 the transformation in 
and out of the rotated frame eliminates the weights of the two 
rates given in Eqs. d35l >. and instead the dynamics is governed 
by the rate: 



ri'fc = tt/2) = 



27V 2 



1 — sine — 



(B7) 



and an additional non-zero contribution from T'2 = 2j. The 
weights of these two rates in the solution for the x component 
are found as: 




(B8) 



where A and B are defined in Eqs. (l34l i. Note that at (p = it/ 2 
there is no dynamics in the y and z axes for our initial state. 

Fig- HH shows the time evolution of the Bloch vector com- 
ponents in the non-rotated frame under 11 -pulse PDD se- 
quence. Fig. [TTT a) captures a weakly-coupled fast TLF with 
parameters as in Figs. |3jb) and (c), while Fig. [TTT b) shows 
a weakly-coupled slow TLF with parameters as in Fig. [3je). 
The dashed red lines depict the transfer matrix calculation 
while the solid blue lines correspond to the analytic solution 
given by Eqs. d27l i, ( |28T i, and dB6l ). In both cases the most 
notable difference with respect to the results of Fig. |3]is the 
prominent dissipative dynamics associated with the z compo- 
nent. The y component (not shown) is an order-of-magnitude 
larger than its counterpart in Fig. [3] but is still considerably 
smaller. 

Similarly, Figs. [TTT c). and (d) show qubit dynamics under 
1 1 -pulse PDD sequence for a strongly-coupled TLF, with pa- 
rameters corresponding to Figs. |4]and|5] respectively. The 
mixing between the x and z components, caused by the rota- 
tion, is apparent in the strong case as well. The y component 
shows appreciable dynamics only when A > v, as seen by the 
dotted red line in Fig. fTTT d). For finite Sh and zero exchange 
(ip = it/2) the x component in all cases shows comparable 
(slightly faster) decay to those given in Figs. [3jf), Efd), and 
[5jc). As seen from Eqs. ( |B6t there is no relaxation in this 
case (y{t) = z(t) = 0). 
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FIG. 11: (color online) Non-rotated Bloch vector components under 
11-pulse PDD sequence for working point 6h = J (tp = 7r /4). (a) 
Weakly-coupled fast TLF (I/7 = 0.1/Lts, « = 2.5 neV), A ~ 7 > 
«; (b) Weakly-coupled slow TLF (I/7 = 0.1 ms, v = 8.4 peV); 
A > 7 > v; (c) Strongly-coupled TLF (I/7 = 0.1 ms, v = 2.5 
neV, Sh = 0.125/xeVO, A > v > 7; (d) Strongly-coupled TLF 
(I/7 = 0.1 ms, v = 2.5 neV, Sh = O.OljuelQ, A > v > 7. 
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